Incommensurate Transverse Anisotropy Induced by Disorder and Spin-Orbit- Vibron 

Coupling in Mn 12 -acetate 
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It has been shown within density-functional theory that in Mni2-acetate there are effects due to 
disorder by solvent molecules and a coupling between vibrational and electronic degrees of freedom. 
We calculate the in-plane principal axes of the second-order anisotropy caused by the second effect 
and compare them with those of the fourth-order anisotropy due to the first effect. We find that the 
two types of the principal axes are not commensurate with each other, which results in a complete 
quenching of the tunnel-splitting oscillation as a function of an applied transverse field. 
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PACS numbers: 75.50.Xx, 71.15.Mb, 75.30.Gw, 75.30.Et 
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The observation of resonant tunneling of magnetiza- 
tion in the single molecule magnet Mni^-acetate 1 * (here- 
after M1112) with a ground-state spin of S — 10 has led to 
many experimental and theoretical investigationsiSi 3 ^ 4 * 5 ^ 
A simple anisotropy Hamiltonian for S = 10 pro- 
vides an excellent approximation to the physics occur- 
ring at low temperatures. Deriving the Hamiltonian 
from density-functional (DF) calculations is challenging 
but possibleA 5 ^ For a particular molecular geometry, a 
magnetic anisotropy tensor may be calculated consider- 
ing spin-orbit coupling within a DF framework. In the 
principal-axes coordinates, the lowest-order spin Hamil- 
tonian can be simplified to the following form: 



Ho = DS 2 z +E(S 2 x -S 2 y ) 



(1) 



where D and E are the uniaxial and second-order trans- 
verse anisotropy parameters and S z is the easy-axis com- 
ponent of the spin operator S. The S4 symmetry of 
the ideal Mni2 molecule causes the value of E to van- 
ish and the lowest-order transverse terms to be fourth 
order. Only transverse anisotropy terms are respon- 
sible for the resonant tunneling between energy levels 
that are almost degenerate. Magnetic tunneling mea- 
surements, however, showed that some of the resonant 
tunneling occurred at a level lower than fourth order* 2 " 
To understand this anomaly in the tunneling, Cornia et 
ali proposed that inherent disorder in solvent molecules 
may break the S4 symmetry and provide nonzero values 
of E. Recent electron paramagnetic resonance (EPR) 
experiments^ and magnetic tunneling measurements^^ 
revealed that the model of Cornia et aZi needed to be 
refined for a quantitative comparison with experiment. 
Our DF calculations^ on the Mni2 showed that the val- 
ues of E for possible configurations of the disordered sol- 



vent molecules are in good quantitative agreement with 
experiment. 

For the Mni2, the calculated D value is -0.556 
which agrees well with experiment. The calculated (mea- 
sured) value of D, however, does not account for all the 
measured anisotropy barrier of 65 One needs a lon- 
gitudinal anisotropy of order higher than second-order. 
In this spirit a fourth-order spin-orbit-vibron (SOV) in- 
teraction was proposed by Pederson et <;iJ-L. Ideally, this 
SOV interaction can only contribute to modifying the 
second-order barrier and to activating two fourth-order 
transverse terms in the spin Hamiltonian. At this level 
a strong SOV interaction can complicate the tunnel- 
ing experiments in several different ways. The simplest 
complication is that there are two different longitudi- 
nal 4th-order terms which scale as S 2 S 2 and S z respec- 
tively. If the S z term is dominant, then at any reso- 
nant field only a single pair of states are involved with 
tunneling. The fourth-order transverse terms result in 
departures of the period of the tunnel-splitting oscilla- 
tion from AH X = 2y/2E(\D\ +E)/g/j,B that has been 
derived by GargiiS. Such departures were first identified 
by Wernsdorfer et alM* and later modified to include 
the 4th-order transverse terms by the group of Garg. 14 
Recent EPR and magnetic experiments^ suggested that 
for Mni2 the in-plane principal axes of the second-order 
anisotropy may not be aligned with those of the fourth- 
order anisotropy so that the oscillation in the tunnel split- 
ting could be completely quenched in this material, in 
contrast to Fe«[l5|. 

In this paper, we perform the DF calculations on a sta- 
tistically weighted collection of three different E values 
caused by the solvent disorder and on the fourth-order 
SOV interaction. We then combine these two effects to 
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determine whether DF theory (DFT) can qualitatively 
account for a suppression of oscillations in the tunnel 
splittings that occur in transverse tunneling experiments. 

Our calculations have been performed within DFT 
using the NRLMOL program^. The NRLMOL calcu- 
lates an accurate Gaussian-type orbital representation 
for the self-consistent occupied and unoccupied molec- 
ular orbitals. We use DFT within the generalized- 
gradient approximation 17 to account for the quantum- 
mechanical behavior of the electrons. The geometries 
of the molecules are fully relaxed. Considering the 
spin-orbit coupling for a relaxed geometry, as discussed 
in Ref. the second-order anisotropy Hamiltonian 
{^a,b— x,y ,zlabS a Sb) may be derived from the calculated 
orbitals. Diagonalization of the 7 matrix provides the 
2nd-order anisotropy barrier and the principal axes. 

As shown in Fig. ^ a single Mni2 molecule is 
surrounded by four acetic-acid solvent molecules 
(CH3CO2H), each of which is shared by two neighboring 
Mni2 molecules. If a Mni2 molecule is hydrogen-bonded 
to four CH3CO2H molecules as pictured in Fig. ^ 
then the symmetry of this molecule remains the same 
as S4. Since CH3CO2H can bind to either of the 
Mni2 neighbors with the same energy 1/16 of the 
Mni2 molecules have the S4-symmetry with hydrogen 
bonds and another 1/16 of the molecules have the S4 
symmetry without hydrogen bonds. The remaining 
molecules have different orientations which break the 
symmetry. As discussed in Ref. there are a total of 
six different configurations which have the number of 
hydrogen bonds n, statistical weights, and D and E val- 
ues of (0,1/16,-0.54,0.000), (4,1/16,-0.56,0.000), 
(1, 4/16, -0.54, 0.008), (2, 4/16, -0.55, 0.000), 

(2,2/16,-0.55,0.016), and (3,4/16,-0.55,0.008). 
The two distinctive cases of two hydrogen bondings 
are characterized by having neighboring acetic acids 
("cis") hydrogen-bonded or next neighboring acetic acids 
("trans") bonded, respectively. The S4 symmetry sug- 
gests that in the case of significantly small off-diagonal 
elements in the 7 matrix for n = 1, the perturbations 
due to the "cis" -geometry should cancel each another 
and that those due to the "trans" -geometry should add 
constructively. Thus, we find that 5/8 of the molecules 
have appreciable E values of 0.008-0.016 K. 

The Hamiltonian for a single uniaxial anisotropic spin 
coupled to a one-dimensional harmonic oscillator is given 
by H = lzz S 2 + \{P 2 + u; 2 Q 2 ) + QEabJab^aSb, where 
P, Q, and lj are the momentum, position, and fre- 
quency of the oscillator and 7 Qb = djat/dQ. As shown 
in Ref. 01 > the energy of the coupled system as a func- 
tion of an eigenvalue of S z (M) is £ = ui/2 + DM 2 — 
{A + BM 2 ) 2 /2uj 2 where A = S(S + l)( 7ra + j' yy )/2 and 
B = 7 22 — (j xx + j yy )/2. Generalizing this problem to 
the Mni2 requires coupling all normal modes. In gen- 
eral, the fourth-order corrections to the Hamiltonian are 
Ti 4 = ?-abcdA abcd S a S b S c S d . Only the Raman active vi- 
brational modes lead to a change in the fourth-order bar- 
rier. 




FIG. 1: Geometry of a M1112 molecule with four hydrogen- 
bonded acetic acids molecules. This high-symmetry configu- 
ration accounts for only 1/16 of the total concentration. For 
details, see Ref.jg. 

Combining the SOV interaction with the symmetry- 
breaking effects of the solvent disorder, we can write the 
total spin Hamiltonian as: 

H = DS 2 + E[cos{2a){S 2 - S 2 y ) + 2sin(2a)S x S y } 

+A 1 S 2 (S 2 - S 2 /3) + A 2 [3S 4 + 35Sj - 3QS 2 S 2 } 
+B 1 (S 4 + S A y - 6S 2 X S 2 ) + B 2 [S x S y {S 2 x - Si)] ,(2) 

where a denotes the angle between the medium axes of 
the calculated second-order and fourth-order anisotropy 
and A\ , A 2 , B\ and B 2 represent fourth-order anisotropy 
parameters. If one replaces the operators (S x , S y , S z ) by 
classical spin vectors and recast the fourth-order trans- 
verse terms in spherical coordinates, then the classical 
potential energy in the xy plane (e.g. 9 = tt/2) becomes 

£ = £ o + S 4 [Si cos(40) + ^ sin(40)] , (3) 

where £q is a constant term and <j) is the azimuthal an- 
gle. This shows that the energy surface of a system 
with S4 symmetry, at fourth order, resembles a four-leaf 
clover in the xy plane. Furthermore, when the fourth- 
order terms in Eq. @ are accounted for, the in-plane 
principal axes are determined to lie along the nodes or 
antinodes of the classical potential. In addition, once 
B\ and B 2 are determined in given coordinates, a ro- 
tation by cos~ 1 [ J Bi/(B 2 + B 2 /16) 1/2 } allows one to re- 
express the two transverse terms as C(S\ + St.) with 
C = [B\ + B|/16] 1 / 2 /2. The experimental value of C is 
2.3 x 10~ 5 K£ while the SOV contributions to B x and B 2 
lead to C sov = 0.053 x lO^Kii 

Even with the rotation of the B\ and B 2 terms to 
the more usual definition, the 2nd-order transverse terms 
may not, in general, be simultaneously expressed in the 
diagonal form. So the dependence on a is required. To 
compute a we calculate the separate magnetic anisotropy 
energies from the SOV interaction and the solvent disor- 
der as a function of <p in the xy plane. For the solvent 
disorder only three configurations [n = 1, 2(trans),and 
3] are considered in the calculation because the rest of 
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FIG. 2: Logarithm of the M = ±5 tunnel splitting versus 
applied transverse field for a Mni2 with one 13 C isotopic defect 
and pure Mni2- Different values for a (0, 3, 12, and 24°) were 
used. 

the three configurations do not produce E. We find that 
a is 24°. This is a bit smaller than the experimentally 
measured value of a = 27° £ We calculate the tunnel 
splitting between M = — 5 and M — 5 with various val- 
ues of a in zero field from diagonalizing the quantum 
mechanical form of Eq. which takes into account the 
non-commutativity of the spin operators. For the pure 
Mni2 sample we find that the oscillation in the tunnel 
splitting is entirely quenched at the calculated value of 
a = 24° as shown in Fig.|^b). When a = 0, as expected 
from Fesr^ the oscillation is clearly visible. Even for a 
very small value of a such as 3°, the amplitude of the os- 
cillation tends to be mostly damped out. As the applied 
transverse field increases, the effect of the incommensu- 
rate transverse anisotropy decreases because the tunnel- 
ing is governed by a strong transverse field. Because of 
the fourth-order terms in Eq. , the period of the oscil- 
lation decreases as the transverse field increases^ 
For a naturally occurring collection of Mni2 molecules, 



37% of the molecules contain one or two isotopic defects. 
In the case of isotopic defects in the Mni2, the SOV in- 
teraction provides additional symmetry-breaking in the 
spin Hamiltonian. If the SOV interaction is large enough, 
this fraction of molecules would contribute to an almost 
continuous distribution of broken-symmetry anisotropy 
Hamiltonians. We calculate the same tunnel splitting 
for the Mni2 doped with one carbon isotope ( 13 C) for 
various values of a [Fig. [3a)] ■ We find that even one 
isotopic defect substantially suppresses the oscillation in 
the tunnel splitting but the period remains unperturbed. 
Since the tunneling experiments were performed in the 
presence of the longitudinal field that makes M and M 1 
almost degenerate where M + M' ^ 0, direct comparison 
with experiment is not possible. The effect of the longi- 
tudinal field and the different types of isotopic defects on 
the tunnel-splitting is in progress. 

In summary, for Mni2 we have considered the second- 
order anisotropy induced by the disordered solvent 
molecules and the fourth-order SOV interaction within 
the DFT framework. The second-order transverse terms 
could cause the oscillation in the tunnel splitting. Exper- 
imental data indicated that the in-plane principal axes 
of the second-order anisotropy was not commensurate 
with those of the fourth-order anisotropy. DFT calcu- 
lations showed that the angle between the two types 
of anisotropy and resulting quenching of the oscillation 
agree well with experiment. 
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